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ABSTRACT 

The MSSM contains CP-violating phases that may have important observable effects in 
Higgs physics. We review recent highlights in Higgs phenomenology obtained with the 
code CPsuperH, a useful tool for studies of the production, mixing and decay of a coupled 
system of the neutral Higgs bosons at future high energy colliders such as the LHC, ILC 
(7LC), and a muon collider (MC). CPsuperH implements the constraints from upper limits 
on electric dipole moments, and may be extended to include other related low-energy 
observables, such as b — > 37 and B — > Kll, and to compute the relic abundance of the 
lightest neutralino. 
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1 Introduction 



CP violation is one of the most sensitive probes for new physics beyond the Standard 
Model, and some additional source of CP violation is needed to account for cosmological 
baryogenesis. Ongoing experiments provide no indications of any flavour or CP-violating 
physics beyond the CKM mixing paradigm. Nevertheless, -B-factories and the Tevatron 
collider are constantly refining their probes of the CKM mechanism, and the LHC and other 
TeV-scale colliders will open up new vistas in the study of CP violation. Supersymmetry is 
one of the more attractive extensions of the Standard Model, and there are many possible 
supplementary sources of CP violation even within the minimal supersymmetric extension 
of the Standard Model (MSSM). These are reduced to just two if one assumes universal soft 
supersymmetry-breaking parameters, which would be sufficient to facilitate baryogenesis 
at the electroweak scale. They would also have a wealth of implications for many different 
areas in particle-physics phenomenology, notably through radiative corrections in the Higgs 
sector [1-7] . CP- violating effects may play a prominent role in Higgs production at future 
high energy colliders, such as the LHC [8-25], ILC [26,27], 7LC [28-36] and MC [37-45]. 
In addition, CP-violating phenomena mediated by Higgs-boson exchanges may manifest 
themselves in a number of low-energy observables such as electric dipole moments [46-61]. 
The CP phases of the MSSM may also affect flavour-changing neutral-current processes 
and CP asymmetries involving K and B mesons [62-85]. Moreover, CP- violating Higgs 
effects may influence the annihilation rates of cosmic relics and hence the abundance of 
dark matter in the Universe [86-89]. Finally, we recall that an accurate determination of 
the Higgs spectrum in the presence of CP violation is crucial for testing the viability of 
electroweak baryogenesis in the MSSM [90-102] 

The Fortran code CPsuperH [103] is a powerful and efficient computational tool for 
understanding quantitatively such phenomenological subjects within the framework of the 
MSSM with explicit CP violation. It calculates the mass spectrum and decay widths of the 
neutral and charged Higgs bosons in the most general MSSM including CP-violating phases. 
In addition, it computes all the couplings of the neutral Higgs bosons #1,2,3 and the charged 
Higgs boson H + . The program is based on the results obtained in Refs. [104-106] and the 
most recent renormalization-group-improved effective-potential approach, which includes 
dominant higher-order logarithmic and threshold corrections, 6-quark Yukawa-coupling re- 
summation effects, and Higgs-boson pole-mass shifts [6, 7] The masses and couplings of the 
charged and neutral Higgs bosons are computed at a similar high-precision level. Even in 
the CP-conserving case, CPsuperH is unique in computing the neutral and charged Higgs- 
boson couplings and masses with equally high levels of precision, and is therefore a useful 
tool for the study of MSSM Higgs phenomenology at present and future colliders. 



2 CPsuperH 

The tarred and gzipped program file CPsuperH. tgz can be downloaded from 



http : / /www . hep . man . ac . uk/ u/ j slee/ CPsuperH . html 

When running CPsuperH, the main numerical output is stored in arrays. These include 
the masses of the three neutral Higgs bosons, labelled in order of increasing mass such 
that Mh 1 < Mjj 2 < Mh 3 - Since the neutral pseudoscalar Higgs bosons mixes with the 
neutral scalars in the presence of CP violation, the charged Higgs-boson pole mass M#± 
is used as an input parameter. CPsuperH also yields the 3x3 Higgs mixing matrix, O a i. 
(0i,0 2 ,a)a = O ai (Hx, H 2 , H 3 )J and all the couplings of the neutral and charged Higgs 
bosons. These include Higgs couplings to leptons, quarks, neutralinos, charginos, stops, 
sbottoms, staus, tau sneutrinos, gluons, photons, and massive vector bosons, as well as 
Higgs-boson self-couplings. We note that, in addition to quantities in the Higgs sector, 
including decay widths and branching fractions, CPsuperH also calculates and stores the 
masses and mixing matrices of the stops, sbottoms, staus, charginos, and neutralinos. For 
a full description, we refer to Ref. [103]. 



3 Collider Signatures 

To analyze CP-violating phenomena in the production, mixing and decay of a coupled 
system [107] of multiple CP-violating MSSM neutral Higgs bosons at colliders, we need a 
"full" 3x3 propagator matrix D(s), given by [108] 

/ s-M^ +z9fmn u (s) i9fmfli2(s) z3mfli 3 (s) \ ^ 

D(s) = s iQmU 21 (s) s-M| 2 +^mn 22 (s) z3mfi 23 (s) 

V z3mn 31 (s) z3mn 32 (s) s - + iQmfl 33 {s) J 

(1) 

where s is the center-of-mass energy squared, Mh 1 2 3 are the one-loop Higgs-boson pole 
masses, and the absorptive parts of the Higgs self-energies Qfrnlljj(s) receive contributions 
from loops of fermions, vector bosons, associated pairs of Higgs and vector bosons, Higgs- 
boson pairs, and sfermions. The computed resummed Higgs-boson propagator matrix has 
been stored in the array DH(3,3). 

The so called tri-mixing scenario has been taken for studying the production, mixing 
and decay of a coupled system of the neutral Higgs bosons at colliders. This scenario is 
characterized by a large value of tan j3 = 50 and the light charged Higgs boson M#± = 155 

Some new features appearing in this write-up, for example, the propagator matrix DH(3,3) and some 
low-energy observables, will be implemented in a forthcoming version of CPsuperH. 



GeV. All the three-Higgs states mix significantly in this scenario in the presence of CP- 
violating mixing. Without CP violation, only two CP-even states mix. For details of the 
scenario, see Refs. [108-111]. 



3.1 LHC 

At the LHC, the matrix element for the process g(\i)g(X 2 ) — > H — > f(a)f(a) can conve- 
niently be represented by the helicity amplitude 

M"(aa;X 1 X 2 ) = ^ 9f ^ S °" (a; \i) g S^S XlX2 , (2) 

where a and b are indices of the SU(3) generators in the adjoint representation and a, a, and 
Ai.2 denote the helicities of fermion, antifermion, and gluons, respectively. The amplitude 
(a; X) g is defined as 

(^;A) fl = E [S!(v r s)+ i XPf(v r s)]D l] (s)(aP f g s Hjff - i g^ ff ), (3) 

i,j=l,2,3 

where 

2 

S!(V§) = E 9f9 S Hj f —F sf (r f )- E Hu.JllJ^r'^O- 
f= b ^ > fj=ti,hMM h 

Pf(V$) = T,9f9 P H J f ^-F pf (r f ), (4) 
f=b,t " L f 



where (3f = Jl — 4m^/s, r x = s/Am 2 x and s-dependent Sf and Pf are scalar and pseu- 
doscalar form factors 2 , respectively. The Higgs-boson couplings to quarks 9 H 'Jif and squarks 
gjj.f*f ., and the explicit forms of the functions -F s /,p/,o are coded in CPsuperH [103]. When 
f = t ,t ,x° , X , e tc ; one can construct CP asymmetries in the longitudinal and/or trans- 
verse polarizations of the final fermions which can be observed at the LHC. The CP asym- 
metries can be defined similarly in the case of other production mechanisms such as 6-quark 
and weak-boson fusions. 

When tan/5 is large, 6-quark fusion is an important mechanism for producing neu- 
tral Higgs bosons. In the centre-of-mass coordinate system for the bb pair, the helicity 
amplitudes for the process 6(A) b(X) — > H — > f(a)f(a) are given by 

M b \aa-XX) = -i^l(a;X) b S a ,S xx , (5) 



where 



(a; X) b = E ( A $> 9H t bb + ^Hfbb) D H (*) (vPf 9 S H] jf ~ WHjff) ■ ( 6 ) 

i,i=l,2,3 



2 These s-dependent gluon-gluon-Higgs couplings are stored in arrays SGLUE(3) and PGLUE(3). 
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Figure 1: The ctoss section o~^ Q ^ [pp(l^lK) — ► t~^~t X] (upper-left panel) at the LHC and its 



associated total CP asymmetry Aq^ = [a 



RR 



LL 



WW 
RR 



+ a 



WW] 
LL \ 



(upper-right panel) as 



functions of &a = $A t = $A b = $A T , where Orr^ll = cr(pp(WW) — >• t^^^^). We have 
considered a tri-mixing scenario with $3 = —10° (dotted lines) and —90° (solid lines). For details 
of the scenario and the CP asymmetry, see [108]. The lower frames are for Higgs bosons produced 
in diffractive collisions at the LHC. The lower-left frame shows the hadron-level cross sections 
when the Higgs bosons decay into b quarks, as functions of the invariant mass M with &a = 90° 
and rapidity y = . The vertical lines indicate the three Higgs-boson pole-mass positions. The 
lower-right frame shows the CP-violating asymmetry when the Higgs bosons decay into r leptons. 
See [110] for details. 

For the tri-mixing scenario, 6-quark fusion is the main production mechanism, and dilutes 
the possibly large CP asymmetry in longitudinally-polarized gluon fusion when f = r. 
This dilution is basically due to the dominance of the 6-quark loop contribution to the 
self energies and the experimental impossibility of distinguishing gluon and 6-quark fusion 
events. In this case, the most promising channel for probing Higgs-sector CP violation 
may be the weak-boson fusion process, which can be separated from gg and bb collisions 
by applying a number of kinematic cuts [112], including the imposition of a veto on any 
hadronic activity between the jets [113-117]: W + W~ — > #1,2,3 — ► t + t~ . The cross section 
a[pp(W + W~) — > H — > t + t~X] lies between 0.2 and 0.6 pb and the CP asymmetry is large 
for a wide range of CP phases, see the two upper frames of Fig. ^ 



Higgs-boson production in an exclusive diffractive collision p + p^-p + Hi+p offers 
unique possibilities for exploring Higgs physics in ways that would be difficult or even im- 
possible in inclusive Higgs production [110]. In spite of the low and theoretically uncertain 
luminosity of the process, what makes diffraction so attractive compared to the inclusive 
processes are the clean environment due to the large rapidity gap and the good Higgs-mass 
resolution of the order of 1 GeV which may be achievable by precise measurements of the 
momenta of the outgoing protons in detectors a long way downstream from the interaction 
point. It may be possible to disentangle nearly-degenerate Higgs bosons by examining the 
production lineshape of the coupled system of neutral Higgs bosons, see lower-left frame 
of Fig. Moreover, the CP-odd polarization asymmetry can be measured when the po- 
larization information of Higgs decay products is available, see the lower-right frame of 

Fig.m 



3.2 ILC 



A future e + e~ linear collider, such as the projected ILC, will have the potential to probe 
the Higgs sector with higher precision than the LHC. At the ILC, the Higgs-boson coupling 
to a pair of vector bosons, gUiW = cpO^i + spO^i, plays a crucial role. There are three 
main processes for producing the neutral Higgs bosons: Higgsstrahlung, WW fusion, and 
pair production, see Fig. |21 




» Hi 





Figure 2: Three main production mechanisms of the neutral Higgs bosons at the ILC: (a) Hig- 
gsstrahlung, (b) WW fusion, and (c) pair production. 



The cross section of each process is given by 

ZHi) = g 2 HtVV a^(M J 



a{e e 
a(e + e~ 

a(e + e~ 



vvHi) 
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lis M 



2 ^.WWf^/r 
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denotes the 



where g Hi mv = sign[det(0)}e ijk g HkVV , v e = -1/4 + sjv > a e = V 4 and a \ 
corresponding production cross section of the SM Higgs boson. As is well known, the WW 
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Figure 3: The differential total cross section da{ ot (e~e + — ► ZHi — ► Zff)/d^/p* multiplied by 
B(Z — » l + l~) = B{Z — ► e + e _ ) + — ► as functions of the invariant mass of the Higgs 

decay products \pj? in units of fb/GeV when f = b (upper-left panel) and f = r (upper-right 
panel). The CP-conserving two-way mixing (P0) and three CP-violating tri-mixing (P1-P3) 
scenarios have been taken. The lower two frames show the CP-violating asymmetries when Higgs 
bosons decay into tau leptons. The asymmetries are defined using the longitudinal (lower-left 
panel) and transverse (lower-right panel) polarizations of r leptons. See [111] for details. 

fusion cross section grows as ln(s) compared to the Higgsstrahlung and becomes dominant 
for large center-of-mass energy yfs. In the decoupling limit, M H ± ^200 GeV, the couplings 
for heavier Higgs bosons gH 2 , 3 vv are suppressed. In this case, for the production of H 2 
and H 3 , the pair production mechanism is active since \gH 2 H 3 v\ — \9mvvl ~ 1- When 
M H ± < 200 GeV, the excellent energy and momentum resolution of electrons and muons 
coming from measurements of the Z boson in Higgsstrahlung may help to resolve a coupled 
system of neutral Higgs bosons by analyzing the production lineshape, see the two upper 
panels of Fig. El 

In the CP-invariant MSSM framework, there is a selection rule that only two CP-even 
Higgs bosons can be produced via the Higgsstrahlung process and only two combinations 
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h+ 



Figure 4: The t + t~ production plane in the Higgs-boson rest frame, in the case when the r 
leptons decay into hadrons h^ and neutrinos. The longitudinal-polarization vector Pl(Pl) 
and the transverse-polarization vector Pt(Pt) with the azimuthal angle a (a) ofr~ (r + ) are 
shown. 



of CP-even and CP-odd Higgs bosons can be produced. In other words, the observation 
of three distinct Higgs bosons in Higgsstrahlung and WW fusion and/or of all three pairs 
of Higgs bosons in pair production could be interpreted as a signal of CP violation in the 
MSSM framework. However, such an interpretation relies on the hypothesis that there 
exist no additional singlet or doublet Higgs fields. To confirm the existence of genuine 
CP violation, one needs to measure other observables such as CP asymmetries. In this 
light, the final fermion spin-spin correlations in Higgs decays into tau leptons, neutralinos, 
charginos, and top quarks need to be investigated. For the Higgsstrahlung process in which 
the produced Higgs bosons decay into a fermion pair //, the differential cross-section for 
the Higgsstrahlung process is given by 



P 



da 
dp 2 
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+ PtPt 



X 



N f g 2 f {3f(p 2 



16tt 2 



{(1 + P L P L )C{{p 2 ) + (P L + P L )Ci(p 2 
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where p 2 and Nf are the invariant mass squared and the color factor of the final-state 
fermions, respectively. The polarization coefficients C{ may conveniently be expressed as 
follows: 
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(9) 



with 



(°) f = £<7iwA;(p 2 ) {<Tpg s H , 



h3 



ff 



(10) 



When f = t and tau leptons subsequently decay into charged hadrons ti* and neutrinos, 



see Fig. HI we have 

Pl = cos9~ , Pt = sin6>~ , a = ; Pl = cos9 + , Pt = sm6 + , a = (p + — it , (11) 

by identifying the polarization analyser for r T as a T = ±h T with h T denoting unit vectors 
parallel to the h T momenta in the r T rest frame. The angles 9^ and ip ± are the polar 
and azimuthal angles of h^, respectively, in the r ± rest frame. This implies that the 
polarization coefficients Cj{p 2 ) can be determined by examining the angular distributions 
of the charged hadrons coming from the r-lepton decays. See the lower panels of Fig. El for 
the CP asymmetries a T L = C%/C{ and a T T = Cl/Cl [111]. 



3.3 7LC 

The two-photon collider option of the ILC, the 7LC, offers unique capabilities for probing 
CP violation in the MSSM Higgs sector, because one may vary the initial-state polariza- 
tions as well as measure the polarizations of some final states in Higgs decays [109]. The 
amplitude contributing to 7(Ai)7(A 2 ) — >• H — > f(a)f(a) is given by 

M H = ^^(^;A 1 )^ CT ^ AlA2 , (12) 
where the reduced amplitude 

to A>& = E mVi) + i\P?(V~s)\ Dy{3) Wf&j, - ig P HjJf ), (13) 

is a quantity given by the Higgs-boson propagator matrix Eq. {IJ combined with the pro- 
duction and decay vertices. The one-loop induced complex couplings of the 77-Hj vertex, 
S2{\fs) and P^{\fs), get dominant contributions from charged particles such as the bottom 
and top quarks, tau leptons, W ± bosons, charginos, third-generation sfermions and charged 
Higgs bosons 3 : 

SHVS) = 2 £ N c Q)g f gl 1Jf ^F sf {r f ) - N c Q)g H ^~F,{r h ) 

f=b,t,xf 2 ' /j=ti,2,6i,2,n,2 h 

V 2 



-9 HiVV Fi{t w ) - g HiH+H _ —^F (r H , 



2M 2 H± 



P?(V§) = 2 £ N c Q 2 f g f g P Hjf ^-F pf (r f ), (14) 

f=b,t,xt: 



. ... mj 



where r x = s/Am 2 , Nq = 3 for (s) quarks and Nq = 1 for staus and charginos, respectively. 
For the explicit forms of Fi and couplings, see [103]. 

One advantage of the 7LC option over e + e~ collisions at the ILC is that one can 
construct CP asymmetries even when Higgs bosons decay into muons and b quarks, by 
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Figure 5: The cross sections (left column) and the CP asymmetries Aq (right column) for the 
processes 77 — > bb (upper) and 77 — ► (lower). The QED continuum contributions to the 

cross sections are also shown. The tri-mixing scenario with $3 = —10° and &a = 90° has been 
considered, making the angle cuts 6^ ut = 280 mrad and 9^ ut = 130 mrad. The three Higgs masses 
are indicated by vertical lines. See [109] for details. 



exploiting the controllable beam polarizations of the colliding photons. In the case of 
identical photon helicities Ai = A2 = A, corresponding to the different combinations of 
helicities of the initial-state photons A and final-state fermions a, we have four cross sections: 



& f _ f f) - Me f^iYyf 



where a , A = ± and 



a;X) f H + 2R(s) 2 F? (a; X) f c + 2R(s) F? 9fe>«<7; X) f H (a; A) 



c 1 



(15) 



(16) 



We have included the QED continuum contribution and an experimental cut on the fermion 
polar angle 6 has been introduced: | cos 6\ < z, and co&eL = z f . For the explicit forms 
3 Thc arrays SPH0(3) and PPH0(3) are used for the s-dcpendcnt 7-7-Higgs couplings. 



in 



of the functions F* f , F 2 Z/ , R(s), and the QED amplitude (a; A)£, we refer to Ref. [109]. In 
this case, one can construct CP asymmetries with no need to determine the helicities of 
the final states as 

Af - ^g=±{Oa+ ~ °a-) ? x 

4a(77 -> //) 

where (7(77 — > ff) = + crL^ + 0+_ + 0"-+) /4. Figure El shows the cross section 

a(77 — > //) and the CP asymmetry A Q when f = b and r. 

By considering the interference of Higgs-mediated resonant and QED continuum am- 
plitudes and the possibly measurable final-state polarizations, one can construct more than 
20 independent observables for each decay mode and half of them are CP odd. This suggests 
that one can investigate all possible spin-spin correlations in the final states such as tau 
leptons, neutralinos, charginos, top quarks, vector bosons, etc., with the goal of complete 
determination of CP- violating Higgs-boson couplings to them. The cases of tau-lepton and 
top-quark final states are demonstrated in [109]. 



3.4 MC 



For the complete determination of the CP-violating Higgs-boson couplings to SM as well 
as supersymmetric particles, a MC is even better than the 7LC. At a muon collider, it 
is possible to control the energy resolution and polarizations of both the muon and the 
anti-muon. Compared to the 7LC case, the center-of-mass frame is known very accurately, 
providing much better resolving power for a nearly-degenerate system of Higgs bosons. 

For example, let us consider the process /i~(A)/i + (A) — > f(a)f(a). When the helicities 
of the muons are equal: A = A, the 7- and Z- mediated processes are suppressed by a 
kinematical factor of the muon mass, and the Higgs-mediated amplitude contributing to 
the process is given by 

M H = -g^gf{vA)l$**S\x, (18) 

where gf — grrif /2Mw- The reduced amplitude is 

3 

where s is the invariant muon-cillider energy squared. The observables obtained by con- 
trolling solely the muon and anti-muon polarizations are given by [41] 



g^gjNcPf 



c! ± C{ 



RR/LL g ^ s 

where the coefficients are 

C ( = l£[ K ^\/l^l/„-\/l* 



16 7T S 



± ) ^15/16 



(20) 



<T=± 



cL 



<r,+)i\' + \{<r,-)l 
~3te E [(<r-)K<r,+) 



a=± 



°16 



^E[K-;+);i 2 -K- 

±9m £ [(<r-)t(<r,+)l 



\/|2 



<T=± 



(21) 



1 I 



390 



,— , 60 

a 

i_, 50 

X 40 
— 

tT* 30 
I 

« 20 

s 

w 10 

-10 



390 



-1 
395 



^-4-r' 
400 



405 



410 

Vs 



395 



400 



405 



410 

Vs 



,— , 2 
ft 1.5 



1 E~ 

0.5 E- 

— 

-0.5 ~ 

-1 r 

-1.5 E~ 
-2 



,— , 0.6 
i— i 0.4 



0.2 - 

— 

-0.2 - 

-0.4 - 
-0.6 




390 395 400 405 410 

Vs 



390 395 400 405 410 

Vs 



Figure 6: The cross sections (cr RR +cr LL )/2 (upper left), cry (upper right), {o RR —cr LL )/2 (lower 
left), and a± (lower right) for the process \i~ —* H —* tt as functions ofy/s. Here, (a RR ±o- LL )/2 
is given in units of fb but the other cross sections are in pb. The solid and dashed lines are for 
<&A = 0° and $>a = — 90°, respectively. 

We note that C( and C{ 6 are CP-odd observables. In Fig. El we show (o~ RR ± cr LL )/2 and 
er when Higgs bosons decay into top quarks / = t, as functions of \fs. The parameters 
are chosen as in Ref. [41], except that tan/3 = 5 and Mh± = 407 GeV. The solid lines are 
for the CP-conserving case with $^ = 0° and the dashed lines are for the CP-violating case 
with $ A = -90°. The masses are M H2 = 400.5 GeV and M Hi = 400.6 GeV in the CP- 
conserving case, but we have a larger mass splitting in the CP-violating case: Mh 2 = 396.2 
GeV and M Ha = 402.4 GeV. The results have been updated, compared to those presented 
in Ref. [41], by including the off-diagonal absorptive parts in the Higgs-boson propagators, 
Higgs-boson pole mass shifts, 6-quark Yukawa-coupling resummation effects, etc. More 
observables could be obtained by considering final-state spin correlations and interference 
effects with the processes mediated by 7 and Z bosons. 
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Figure 7: The Thallium EDM dxi = ^ti x 10 24 e cm in the tri-mixing scenario. The upper-left 
frame displays \dr\\ in the (<&a, $3) plane. The unshaded region around the point $3 = &a = 180° 
is not theoretically allowed. The different shaded regions correspond to different ranges of \dr\\, 
as shown: specifically, |rfri| < 1 in the narrow region denoted by filled black squares. In the upper- 
right frame, we show |g!ti| as a function of &a for several values of $3. In the lower-left frame, 
we show \d^i\ as a function of $3 for four values of &a- In the lower-right frame, we show the 
Cs (dotted line) and d e (dash-dotted line) contributions to df \ separately as functions of $3 when 
= 60°. As shown by the dashed line, the chargino contribution is negligible. See [111] for 
details. 



4 Low-energy Observables 



Low-energy observables provide indirect constraints on the soft SUSY breaking parameters. 
The observables are particularly useful for identifying the favoured range of parameter space 
when the SM predictions for them are strongly suppressed and/or precise experimental 
measurements of them have been performed. Such observables include EDMs, (g — 2)^, 
BR(b -> s 7 ), A C p(b -> 57), BR(B -> Kl I) and BR(B s d -> 

Currently, the EDM of the Thallium atom provides one of the best constraints on 
the CP-violating phases, depending on the SUSY scale. The main contributions to the 
atomic EDM of 235 T1 come from two terms. One of them is the electron EDM d e and the 
other is the coefficient Cs of a CP-odd electron-nucleon interaction. The coefficient Cs 
is essentially given by the gluon-gluon-Higgs couplings and the two-loop Higgs-mediated 



13 



electron EDM [58,59] is given by the sum of contributions from third-generation quarks and 
squarks and charginos 4 . We show in [111] that it is straightforward to obtain C$ and the 
Higgs-mediated d e by use of couplings calculated by CPsuperH. As can be seen for example 
in Fig. one can then implement the thallium EDM constraint on the CP- violating phases 
and demonstrate that they leave open the possibility of large CP-violating effects in Higgs 
production at the ILC. 

The code CPsuperH will be extended in near future to include CP-violating effective 
FCNC Higgs-boson interactions to up- and down-type quarks [85,118-120]. The deter- 
mination of these effective interactions may be further improved in the framework of an 
effective potential approach, where the most significant subleading contributions to the 
couplings can be consistently incorporated. At large values of tan/5, Higgs-mediated inter- 
actions contribute significantly to the £>-meson observables mentioned above and so may 
offer novel constraints on the parameter space of constrained versions of the MSSM, such 
as the scenario of minimal flavour violation. 

5 Cosmological Connections and Future Directions 

Imposing R parity, the Lightest Supersymmetric Particle (LSP) is stable and becomes one 
of the strongest candidates for the cold dark matter which provides some 23 % of the 
energy density of the Universe. When the lightest neutralino is the LSP, the annihila- 
tion of neutralino pairs and scattering off ordinary matter become relevant in Cosmology. 
Neutralino-pair annihilation cross sections are needed in calculations of the cosmological 
neutralino relic density and of the fluxes of cosmic-ray antiprotons, positrons and gamma- 
rays due to relic annihilations in the Galactic halo. The fluxes of high-energy neutrinos 
from the cores of the Sun and Earth depend on the scattering processes leading to neu- 
tralino capture, as well as on the rates for different annihilation channels. Scattering cross 
sections are also important for the direct detection of the neutralino via its interaction with 
matter. 

In the presence of CP-violating mixing in Higgs sector, the couplings of Higgs bosons 
to the SM and SUSY particles are significantly modified and the contributions to the 
pair-annihilation and scattering cross sections from Higgs mediated processes could be 
largely affected. For example, for non-relativistic neutralinos, the amplitude of the process 
XiXi "~ Hi — > WW , H\H\ is highly suppressed in CP-invariant theories, but may be 
enhanced when the intermediate Higgs bosons are CP-mixed states [89]. We note that 

4 Hcrc we have not considered one-loop contributions, since they are independent in the absence of 
any assumptions relating the soft SUSY-breaking terms of the first two generations to those of the third 
generation. 
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CPsuperH has been used for the calculation of the relic abundance of the lightest neutralino 
in Ref. [121]. 

As mentioned above, CP- violating phases in the MSSM may also make possible baryo- 
genesis at the electroweak scale [100-102, 122]. A key ingredient in verifying this and other 
cosmological connections will be laboratory measurements of CP- violating phenomena, and 
CPsuperH will be a key tool for evaluating their implications for the underlying model pa- 
rameters by extending its scope to low-energy and cosmological observables. 
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